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Abstract. 

The phase transitions in the Bose-Hubbard model are investigated. A single-particle Green's func- 
tion is calculated in the random phase approximation and the formalism of the Hubbard operators is 
used. The regions of existence of the superfluid and Mott insulator phases are established and the (fJL,t) 
(the chemical potential - transfer parameter) phase diagrams are built. The influence of temperature 
change on this transition is analyzed and the phase diagram in the (T,fi) plane is constructed. The 
role of thermal activation of the ion hopping is investigated by taking into account the temperature 
dependence of the transfer parameter. The reconstruction of the Mott-insulator lobes due to this effect 
is analyzed. 
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INTRODUCTION 

The Bose-Hubbard model (BHM) has been intensively investigated in the last 15 years. This model 
is of great interest due to the experimental realization of Bose-Einstein condensation of ultra-cold 
atoms in optical lattices (see, for example, yj]). Various theoretical methods were used to study 
the model: mean-field theory [2], random phase approximation lElH], strong coupling approach 
A, quantum Monte-Carlo method |@], bosonic version of dynamical mean field theory 0]. The 
Hamiltonian of the model includes two terms. One is connected with the on-site energy U, which 
describes repulsion of bosons at a lattice site; another one describes the nearest-neighbor hopping 
with tunneling parameter t: 

H = -£^+£ J + y£n^-l)-ML«n (1) 

ij i i 

\i is the chemical potential. The competition between the kinetic and on-site repulsion terms defines 
the equilibrium state of the system. The existence of the superfluid and Mott-insulator phases is a 
characteristic feature of this model. When the kinetic energy dominates (t/U ^> 1) the ground state 
of the system is superfluid (SF). In the opposite case the ground state is a Mott insulator (MI). 

In the limit U — > °° this model reduces to the hard-core boson model. Models of a lattice gas 
type with Pauli statistics of particles are often used for the description of ionic conductors and 
calculation of their conductivity starting from the works of Mahan [8] and others Is, 10J. Recently 



the one-particle spectrum has been investigated in the one-dimensional limit fl 1 Qtl . As noted above, 
the Bose-Hubbard model can be directly applied to the description of optical lattices. Besides, the 
model of Bose-Hubbard type can be also usefull not only for investigation of ionic conductivity in 



crystalline ionic conductors, but it can also be applied to describe the intercalation in crystals 
and kinetics of ionic adsorption on the crystal surfaces lll2n . 

In this work we consider the case of the finite value of the on-site interaction U and study 
phase transitions in the BHM at finite temperatures. The most previous investigations of the Bose- 
Hubbard model were restricted to the case of zero temperature. Recently only the low-temperature 
case (T <C U) has been investigated |0,E3l. A special attention will be paid here to the investigation 
of the influence of thermal activation of the ion hopping on the shape of phase diagrams (we take 
into account the temperature dependence of the transfer parameter). The effect of thermal activation 
can be important in the case of ionic conductors. It was not studied in the framework of the BHM; 



such an investigation has been performed up to now only in the hard-core boson limit IllOh . 

THE RANDOM PHASE APPROXIMATION 

We introduce the on-site basis \n >,■ (n is a number of bosons on a site i) and use the Hubbard 
operators X" m = \n >,< m| ; . The Bose-operators of creation and annihilation can be expressed in 
terms of the X- operators: 

^ = £v^TTxf'" +1 , n / = £nXr. (2) 

n n 

Then the Hamiltonian (OQ) can be written as follows: 

H = Hq+Hi, 

Hq = £A„X™; Xn = ^n{n-l)-nn, 

in ~ 

Hi = -J^Ujb+bj. (3) 

ij 

To find the Green's function ((b\b + )) we employ the equation of motion technique: 

co((X^ m+l \X^)) = ^5 Zp ^(xr-X; +1 ' m+1 ) + (^ + i-A m )((X z m ' ra+1 |x; +1 ''-)) 

ij 

- 5^5 /y ((^(V^^f ! - m+2 -V^xf'- 1 '' n+1 )|x; +1 0). (4) 

ij 

We use the random phase approximation (which is an analogy to the Hubbard-I approximation 
for the case of the fermionic Hubbard model) and perform the following decoupling: 

( ( (x™> m - x™ + 1 ' ra+ 1 ) bj |x; +1 ■'') ) » (x mm - x m+ 1 -< m+ 1 ) ( (bj \x r + Lr )) 

((b+X'^ m+2 \X r p +1 ' r )) w 0. (5) 

It should be noted that here we consider the Mott-insulator phase, where (&,•) = 0. If we want to 
investigate the superfluid phase we should take into account that in this phase (bj) ^ (this case is 
considered below). 



As a result, the following expression for the Green's function is obtained: 

A m = — ^/7j, (X mm —X m+ ,m+ ) = Q m . (7) 

We pass to ^-representation and the Green's function G^(o) = ((b\b + ))i cco can be written as 

1 *° 

G k (co) 



We should note that in the independent subband approximation (this approximation is valid in the 
case U ^> t and in the regions of the integer value of fx/U) the Green's function takes the form 

G = 1 y (m+\)Q m 
h 2n^co-A m +t k (m + \)Q m ' 

and the subband energy is equal to X m = A m — (m + l)Q m tk- 

In order to study the phase with the Bose-Einstein condensate we can also apply decoupling of 
the © type, which in this case is more complicated: 

<<(^ / ^-js s : / ^+ 1 ^+ 1 )z, y -|js s :-+ 1 ^>> ^ ^<<z ?y -|^+ 1 ^>> <z, y -><<(^r'--^+ 1 ^+ 1 )|^+ 1 ^>>. (io) 

The function (((X m,OT — X w+ ' m+ )bj\X p + ,r )) satisfies the following equation: 

co((X™> m -X™ +1 > m+1 \X r p +1 > r )) = ±S^(X m > m - l }8 m)r+l + (X m+2 > m+l }8 m+hr 
-2(X-+ 1 '-)6 rar )-£^( v ^((^ m - 1 ^-^r 1 ' m |^ +1 ' r )^ 
+V^+2((X', n+2jn+1 bj - bix™ +1 > m+2 \x r p +1 > r )) 

-2 V / ^TT( (x™ +1 > m bj - bjx'^ m+l . (ii) 

Using decoupling (flOl) in the independent subband approximation (which is equivalent to taking 
into consideration the transition between levels \m) and \m+ 1) only) the equation takes the more 
simple form 

CD((xr i -X^ +hm+l \X r p ^ r )) = --5 lp 5 mr ((X'" +1 n 

+2(m+ l)^t,j((Xp m+l )((X^ +l ' m \X r p + ^ r )) - (X™> m+1 )((Xj l+1 > m \X r p +1 ' r ))) 
j 

-2(m+ \)Y,tij({Xf +l ' n ){{X^ m+l \X r p ^ r )) - (X? +l ' m ){(Xj l ' m+l \X r p +1 > r ))). (12) 



In a similar way, we can obtain equations for the Green's functions ((X w+ ' m \X p +1,r )) and 

+A m )((x^ m \x r p +Ur )) = (m+i)£^((x; +1 ' m |x; +1 -)) 



{{x^ m+1 \x r p +1 ' r )y. 



CO 



J 

+( m +l)^t lj (X l f +l ' m )((X^ m -X^ +l ' m+l \X p ^ r )). (13) 



(co-A m )((X^ m+1 \x;+ l - r )) = ±-8 lp 8 mr Q m - (m+ l^Q^X?^ 1 ^ 1 ')) 

-(m+ l)£? /j (X; J ' m+1 )(((X / ra ' m -X / ra+1 ' ra+1 )|X;+ 1 ''")). (14) 
j 

The levels \m >, \m+ 1 > form two-state subspace where we can introduce pseudospin formal- 
ism: 

v m,m+l _ c — . v m+l,m _ c + , 1 / v m+l,m+l v mm\ _ C Z 

A / - b l,m' A l l,m' 2 ' ' ~ ' ^ 

The mean value <2„, = — 2(5^) is connected with the longitudinal component of the pseudospin, 
while the order parameter (bj) in the condensate phase is proportional to the transverse component 

(bj) « v^+T(5-) = v^+T(sy. 

The set of equations (fT2l) . (TT3T ), and (fl4l) can be rewritten using pseudospin Green's functions. 
The equations are analogous to those obtained in the RPA for the XXZ model. Their solution can 
be written in the form given in nl4\ 

- 1 ^ z \ E m(cos 2 e m + 1) + 2 co cos e m -4(o^}J qyfn cos 2 m ^ 

2,K CO 2 {Em ^{^m)Jq,m){E m 1(o m )Jq m )cOS 2 Q m 

In the mean field approximation the angle 6 m and parameter (o^) are defined by the following 
equations: 

A m sin 6 m - 2(m + l)t (S x m ) cos 6 m = (17) 
1 U PE„ 



(a m ) = 2 tanh ^; E ™ = V A - + 4(m+l)% 2 (5-)2; 
to = t q=0 ; J q>m = (m+l)t q ; (S m ) = -(o£) sin0 ra . (18) 

The solution sin0 m = ((Sf n ) = 0) corresponds to the normal (nonsuperfluid) phase; in the phase 
with condensate cos G m = ~ 2 (m+fj" ((r ) an( * ( G m) * s gi ven by the solution of the following equation 

(^) = itanh(j8(m+l> (^». (19) 
The transition to the phase with condensate takes place at temperature 

I _|_ Asl _j 

K = A m (\n—^) . (20) 



m m \ l_Am 

Jo,m 



This temperature reaches its maximum value T£ ax = Jo, m /2 at A m = mU — \l = 0. 

The case of pair of the states |0 >, 1 1 > corresponds to the hard core boson model; the peculiar- 
ities of the spectrum and density of states were investigated in [11411 . The thermodynamics of the 
model was also studied in jl4 \ in the case of the thermal activated hopping. 



RESULTS 



The phase transition from the Mott-insulator to superfluid phase is characterized by divergence 
of the Green's function G^ =0 (co = 0) — > °°. This condition can be obtained using the following 



reasoning. Let us find the expression for the commutator [X™ ,m+ ,H] in the Hubbard-I type 
approximation: 

[X^ m+ \H] « A m X;' hm+l - Y,tijV^TT(X mm -X m+1 ' m+1 )bj. (21) 

The mean value of this commutator should be equal to zero (because the average (x m ' m+1 ) does 
not depend on time in a stationary state). We get after decoulping: 

TYl ~\~ 1 

<**>[!-**£ -7— Qm]=0. (22) 



in 



Thus, we can see that the nonzero solution (b^) ^ appears in the point of divergence of the 
Green's function Gfc(ct) = 0) — > °° ©. The condition (|22l) in the case £ > (considered in this 
paper) is fulfilled at the center of the Brillouin zone when k = (an uniform phase case) and 
G^ = o(tt) = 0) — ► °° in the phase transition point. In the case f < this condition is satisfied at the 
Brillouin zone edge when k — (71,71,...) (a modulation phase case). 
In our numerical calculations we used the diagonal Hamiltonian Hq to calculate the averages 



(X mm ). In this case the condition (1221) coincides with that obtained in IU3N in the framework of 
the diagrammatic hopping expansion at the calculation of temperature Green's functions (the first 
order in the self-energy in the notations used in Ill3ll ). 
Let us consider the influence of thermal activation. We rewrite the parameter of the ion hopping 



in the following form: t = to&xp(— /3A). Such an approximation was used in [allOfl to take into 
account the presence (in ionic conductors) of a static barrier over which a particle must hop. It can 
be obtained also as a result of renormalization due to the interaction with phonons. 




FIGURE 1. The (W,fl) phase diagram at A = (a) and A ^ (b). 

In Fig. [T] the (W, /i) phase diagram is shown (W is a half width of the initial energy band, 
— W < t < W). The system can pass from the MI to SF phase at the change of the chemical 
potential as well as at the change of the transfer parameter. The influence of thermal activation 
on phase transitions is illustrated in FigflJ?. Thermal activation leads to expanding of the region of 
the Mott-insulator phase. 

The (r,ju) phase diagrams are given in Figf2l The influence of thermal activation is illustrated 
in Figf2b in the case of the one SF lobe Qi/U « 1). At the increase of the activation energy this 




region narrows and then the superfluid phase disappears (there is a critical value A cr ). It should 
be noted that at the fixed value of the chemical potential there are two critical temperatures 
T Cl and T C2 at which the phase transitions between the MI and superfluid phase occur (it takes 
place at A 7^ 0). One can think that the lower critical temperature T C2 can be connected with the 
temperature of the transition to the superionic state in the case of superionic crystals. The simple 
estimate, based on the typical values of t and A parameters for proton conductors (A ~ 0.1... 0.4 
eV; t ~ 0.2... 0.5 eV) shows that T C2 can be of the order of several hundreds of Kelvins, that 
corresponds to temperatures observed in the known family M3H(X04)2 (M=NH4, Rb, Cs; X=S,Se) 
of hydrogen-bonded superionic crystals | |15|1 . 




FIGURE 3. The (U,p) phase diagram at A = (a) and A ^ 0, T = 0.5W (b). 

From the (£/, jii) phase diagrams shown in FigfJ] we can see that the influence of thermal 
activation on the phase diargams is similar to the influence of temperature change in the case 
A = 0. 

In Figj4]the (T,[l) phase diagram calculated in the independent subband approximation using 
the equation (|20l) (where 7o,m = (m+ l)?o ex P( — P^"A)) is shown (dashed line). Solid line is the line 



shown in Figj2^ in the case W = 0.07U. We can see that the independent subband approximation is 
applicable at small values of m. In this region the expression (|20l) describes adequately the nearly 
linear growth of T™ with the increase of the number m. 
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FIGURE 4. Dashed line denotes the phase transition line obtained in the independent subband approx- 
imation using (l20l ). Solid line denotes the phase transition line obtained using (l22l ). W = 0.07£/,A = 0. 



CONCLUSIONS 

The phase transitions in the Bose-Hubbard model at finite temperature are investigated. Models of 
such type can be applied for the description of optical lattices, ionic conductivity, intercalation in 
crystals, and kinetics of ionic adsorption on the crystal surfaces. 

A single-particle Green's function is calculated in the random phase approximation (which is 
an analogy to the Hubbard-I approximation in the case of the fermionic Hubbard model) and 
the formalism of the Hubbard operators is used. This Green's function is obtained for the Mott- 
insulator and superfluid phases (in the last case we use the independent subband approximation, 
which is appropriate for the case of the noninteger (intermediate) mean number of the particle 
concentration per lattice site). 

The regions of existence of the superfluid and Mott insulator phases are established and the phase 
diagrams in the plane (/!,?) (the chemical potential - transfer parameter) are built. The influence 
of temperature change on this transition is analysed and the phase diagram in the (T,ju) plane is 
constructed. The role of thermal activation of the ion hopping is investigated by taking into account 
the temperature dependence of the transfer parameter (t = ?oexp(— /3A)). It is shown that thermal 
activation leads to the reconstruction of the Mott-insulator lobes (in particular, to the narrowing 
of the SF phase region). It is revealed that due to thermal activation of the ion hopping, there are 
two critical temperatures at which the MI-SF phase transition takes place (at the fixed value of the 
chemical potential). A simple estimation shows that the lower temperature can be quite high in 
proton conductors and can correspond to the temperature of the superionic phase transition in such 
compounds. 



REFERENCES 



1. M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, I. Bloch, Nature 415, 39-44 (2002). 

2. K. Sheshadri, H. R. Krishnamurthy, R. Pandit, T. V. Ramakrishnan, Europhys. Lett. 22, 257-263 
(1993). 

3. S. Konabe, T. Nikuni, M. Nakamura, Phys. Rev. A 73 , 033621-12 (2006). 

4. Y. Ohashi, M. Kitauri, H. Matsumoto, Phys. Rev. A 73, 033617-5 (2006). 

5. J. K. Freericks, H. Monien, Europhys. Lett. 26, 545-550 (1994). 

6. B. Capogrosso-Sansone, S. G. Soyler, N. Prokof 'ev, B. Svistunov, Phys. Rev. A 77, 015602-4 (2008). 

7. K. Byczuk, D. Vollhardt, Phys. Rev. B 77, 235 106-14 (2008). 

8. G. D. Mahan, Phys. Rev. B 14, 780-793 (1976). 

9. T. Tomoyose, /. Phys. Soc. Jpn. 66, 2383-2385 (1997). 

10. I. V. Stasyuk, I. R. Dulepa, Condens. Matter Phys. 10, 259-268 (2007). 

11. T. S. Mysakovych, V. O. Krasnov, I. V. Stasyuk, Condens. Matter Phys. 11, 663-667 (2008). 

12. P. D. Reilly, R. A. Harris, K. B. Whaley, J. Chem. Phys. 95, 8599-8615 (1991). 

13. M. Ohliger, A. Pelster, cond-mat.stat-mech/081 0.4399, 
URL [http: //arxiv.org/pdf /0810 . 43 9 9| 

14. I. V. Stasyuk, I. R. Dulepa, Journ. Phys. Studies 13 (2009), in print. 

15. I. V. Stasyuk, N. Pavlenko, J. Phys.: Condens. Matter 10, 7079-7090 (1998). 

Qa30Bi /4arpaMH MO^ejii Bo3e-Xa66ap/j,a npn ciriirceHHiH 

TeMnepaTypi 

LB. CTacioK , T.C. MiicaKOBHM 

iHcmumym $i3UKU KondeHcoeanux cucmeM, 79011 JIbeie, eyA. CeeHVfiupKozo, 1 

AHOTan,ia. ^ocjiiflHceHO cpa30Bi nepexoflii b MO^ejii Bo3e-Xa66ap,a;a. OflHonacTHHKOBa (pyHKirm 
TpiHa po3paxoBaHa b Ha6jiHHceHHi xaoTiiHHiix cpa3 3 BHKopiiCTaHHHM cpopMajii3My onepaTopiB 
Xa66ap/];a. BcTaHOBjierro o6jiacTi icHyBaHHH Ha/rnjiiiHHoi cpa3ii Ta (pa3ii MoTTiBCbKoro 
fliejieKTpHKa i no6y,i];oBaHO (pa30Bi fliarpaMii b njioiuFffli (n,t). IlpoaHajii30BaHO BnjniB 
3mihh TeMnepaTypn Ha dpa30Bi nepexo^n Ta no6y,i];oBaHO dpa30Bi ,a;iarpaMii b njionniHi 
(T,jj,). ^ocjiiflJKeHO BinaflOK TepMinHoi aKTHBar^ii nepecKOKy nacTiiHOK Ta npoaHajii30BaHO 
nepe6y,i];0By cpa30Biix fliarpaM y ubOMy BimaflKy. 

KjiMHOBi cjiosa: MO^ejib Bo3e-Xa66ap,n;a, cpa30Bi nepexo,a;H 
PACS: 03.75.Lm, 05.30.Jp 



